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^SJ ■ Abstract. We present near-IR (J and Ks) number counts and colors of galaxies detected in deep VLT-ISAAC 

' images centered on the Chandra Deep Field and Hubble Deep Field-South for a total area of 13.6 arcmin^. The 

I limiting surface brightness obtained is Ks~22.8 mag/arcsec'^ and J~24.5 (Icr) on both fields. A dlogN/dm relation 

. with a slope of ~ 0.34 in J and ~ 0.28 in Ks is found in both fields with no evidence of decline near the magnitude 

' limit. The median J-Ks color of galaxies becomes bluer at magnitudes fainter than Ks~ 18, in agreement with 

fH , the different number counts slope observed in the two bands. We find a fraction (< 5% of the total sample) 

^ of sources with color redder than J-Ks=2.3 at magnitudes Ks> 20. Most of them appear as isolated sources, 

possibly elliptical or dusty starburst galaxies at redshift z > 2. The comparison of the observed number counts 
with models shows that our J-band and Ks-band counts are consistent with the prediction of a model based on 
a small amount of merging in a f2 = 1 cosmology. On the other hand, we fail to reproduce the observed counts if 
I we do not consider merging independently of the parameters defining the universe. 

^ ' Key words, cosmology - galaxies - 

5t 1- Introduction vations were difficult to reconcile without number evolu- 

. . , . .„ tion (Broadhurst et al. 1992; Gronwall & Koo 1995; Babul 

Near infrared (IR) observations provide significant advan- ^ Rees 1992) 

tages over optical ones in studying galaxy evolution be- 
cause of the small and almost galaxy type independent Disentangling the effects of luminosity and number 
k-corrections at these wavelengths and of their insensitiv- density evolution is fundamental in understanding the na- 
ity to star formation activity over a wide redshift range ture of field galaxy formation and evolution allowing a 
(e.g. Cowie et al. 1994). They also provide a first guess direct comparison with models. For instance, hierarchi- 
about the stellar mass content of the galaxies since they cal models, which assume that galaxies assemble through 
trace the underlying old stellar population of galaxies out merging from smaller sub-units, make concrete predic- 
tor 2-3. These are the main reasons why deep near-IR tions about the evolution of the merger rate, morpholog- 
selected samples have been thought to provide constraints ical mix and redshift distribution (e.g. Baugh et al. 1996; 
on the different galaxy formation and evolution scenarios. Kauffmann 1996; Kauffmann & Chariot 1998). In this see- 
When the first large IR detectors became available, nario the number density of massive (and luminous in the 
near-IR galaxy counts were initially considered a power- near-IR) galaxies should decrease with increasing redshift. 
ful cosmological test (Cowie et al. 1990; Gardner et al. Consequently, an IR selected sample is better suited to in- 
1993). It was later realized that optical and near-IR obser- vestigate the density of these galaxies. 

Send offpnnt requests to: Paolo Saracco, e-mail Kauffmann & Chariot (1998) suggest that the redshift 

saracco@merate.mi.astro.it distribution of a K-band selected sample can provide an 

* Based on observations collected at the ESO-VLT telescope observational test able to discriminate between hierarchi- 

(Prog. ID.64.O-0643, PI E. Giallongo; Prog. ID.164.O-0612, PI cal models and pure luminosity evolution (PLE) models. 

M. Franx). They show that hierarchical models predict a substan- 
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tially lower fraction of high-z galaxies at any redshift than 
do PLE models. Moreover, PLE models predict a high 
redshift tail of the redshift distribution of galaxies much 
more extended than in hierarchical models. Fontana et al. 

(1999) applied this test to a composite sample, deriving 
the photometric redshift distribution of galaxies down to 
K< 21. They found only ~ 5% of galaxies at z > 2 in 
good agreement with predictions of hierarchical models. 

Near-IR colors have been used, e.g. by Eisenhardt et al. 

(2000) , to select z > 1 galaxies with the aim to constrain 
the fraction of high redshift galaxies in a K-band selected 
sample. 

We have obtained deep J and Ks band VLT-ISAAC 
observations centered on the Chandra Deep Field (~ 6 
arcmin^, CDF hereafter) and Hubble Deep Field South 
(~ 7.4 arcmin^, HDFS hereafter) which extend to J'^ 24.5 
and Ks^ 23. The HDFS data are in common with those 
of a similar proposal. This data set constitutes the widest 
area surveyed at these near-IR magnitude limits. In this 
paper we present counts and near-IR colors of galaxies 
detected in these two fields. The plan of the paper is as 
follows: in §2 we describe the observations and the data re- 
duction and show in §3 the procedure adopted to construct 
our photometric samples and to derive number counts. In 
§4 we present color distributions and in §5 we compare the 
observed number counts with different models. In §6 we 
use color selection to identify z > 1 and z > 2 galaxies and 
to obtain lower limits to their number densities. Finally, 
in §7 we summarize our results and conclusions. 

2. Observations and Data Reduction 

2.1. Observations and Photometric Calibration 

The data have been obtained with the ISAAC infrared 
imager/spectrometer (Moorwood et al. 1999) at the ESO 
VLT-UTl telescope. ISAAC is equipped with a 1024 x 1024 
pixel Rockwell Hawaii array providing a plate scale of 
0.147 arcsec/pix and a total field of view of about 
2.5x2.5 arcmin. The two observed fields are centered at 
03h32ml6s, -27°47'25" the center of the Chandra Deep 
Field and at 22h32m55s, -60°33',08" the center of the 
Hubble Deep Field-South. The observations were gathered 
over several nights from September to December 1999 un- 
der homogeneous seeing conditions: around 0.6 arcsec in 
the case of the HDFS and around 0.7 arcsec in the case 
of the CDF. The CDF was imaged in the J and Ks bands 
while the HDFS in the Js, H and Ks bands (see Tab.l). 

The standard "auto-jitter" mode was used to take all 
the images, with the telescope being offset by random 
amounts up to 30 arcsec (CDF) and 20 arcsec (HDFS) 
between individual short exposures. In Tab. 1 we summa- 
rize the number of frames and the relevant total exposure 
times in each band for the two observed fields. 

Photometric calibration of the observations has been 
made by observing several standard stars from the list 
of Infrared NICMOS Standard Stars (Persson et al. 1998) 
with magnitudes ranging between 10 and 12. Instrumental 



total flux has been estimated by deriving the growth curve 
for each star using the IRAF task phot. The estimated 
magnitudes have then been corrected for atmospheric ex- 
tinction assuming Aj=0.1, Ah=0.04 and Aks=0.05. The 
typical uncertainty in the derived photometric calibration 
ranges between 0.02 and 0.04 mag. 

2.2. Data reduction 

Raw frames have been first corrected for bias and dark 
current pattern by subtracting a median dark frame. The 
fiat field corrections have been made by using a mean dif- 
ferential sky- light fiat obtained by averaging the difference 
between a set of low-count sky flat images from a set of 
high-count sky flat images. 

Frames have been then processed by DIMSUM|^ to pro- 
duce the sky subtracted frames on the basis of the follow- 
ing recipe. A first pass of DIMSUM produces a "raw" 
final co-added image by which to derive a master mask 
frame flagging pixels belonging to sources. This master 
mask is de-registered to create a proper mask relevant to 
each frame. For each frame a sky background image is thus 
generated by averaging a set (from 3 to 6) of time adjacent 
frames where sources are now masked out. This procedure 
makes it possible to obtain sky subtracted frames where 
sources are not surrounded by dark halos typical of a sky 
subtraction process performed without masking objects. 

The sky-subtracted frames have been then inspected 
and sky residuals were removed where present. Sky resid- 
uals have been modeled using the IRAF task imsurfit by 
fitting a bi-cubic spline to the frame after masking out the 
sources. 

Frames have been then rescaled to the same airmass 
using the airmass information stored in the header of 
the images and to the same zeropoint by comparing the 
photometry of all the sources detectable in each frame. 
Finally, before shifting and co-adding, frames have been 
additively rescaled to the same median value. 

Shifting and co-adding have been then performed with 
DIMSUM. The final co-added image is the average of the 
shifted frames. The quality of the final images is shown in 
Table 1 where the measured FWHM and the la limiting 
surface brightness are reported. 

3. Number Counts 

3.1. Object catalogs and magnitudes 

Object detection has been performed using the SExtractor 
package (Bertin and Arnouts 1996). We used a Gaussian 
function with a FWHM matching the one measured on 
the frames to convolve the image and an RMS weighting 
map to optimize the detection procedure. Since the two 
data sets have been obtained under slightly different see- 
ing condition, we used different detection thresholds such 

^ Deep Infrared Mosaicing Software, a package of IRAF 
scripts by Eisenhardt, Dickinson, Stanford an d Ward, avail- 



able at 'tp://iraf.noao.edu/contrib /dimsum V2| . 
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0.16 
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Ks 


2.16 


0.27 
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29220 
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Table 1. Details of observations and image quality 
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Fig. 1. The figure compares the Ks-band counts obtained in this work in the CDF and in the HDFS with those in 
the literature. The Maihara et al. (2000) counts have been plotted down to their S/N^^S limit. The slope at Ks> 19 
is 0.28 both for the CDF and the HDFS Ks-band counts. 



that the faintest detectable source in both the fields has 
a minimum signal-to- noise ratio of S/N=3 over the seeing 
disk. The raw catalogs were then been cleaned of those 
objects having saturated pixels and/or lying too close to 
the edges to perform a reliable flux estimate. 

The Ks-band selected catalogs thus obtained con- 
tain 332 sources over an area of 6.01 arcmin^ (CDF) 
and 414 sources over 7.47 arcmin^ (HDFS). The 
CDF near-IR catalog is available in electronic form 



via http: / /www. merate.mi.astro.it / ^saracco . The multi- 
band HDFS catalog will be published in a forthcoming 
paper (Vanzella et al. 2001). 

In order to establish the most reliable and unbiased 
estimate of the "total" flux of sources in our images, we 
compared different magnitude estimates among them. We 
focused our attention on three different issues: 1) the mag- 
nitude which minimizes the systematic underestimate of 
the flux at faint magnitudes; 2) identifying the most re- 
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Fig. 2. The J-band galaxy counts obtained in this work are compared with those in the literature. The Maihara et al. 
(2000) counts have been plotted down to their S/N^3 limit. Our counts continue to rise with a power law slope at 
J> 20 of 0.37 and 0.34 in the CDF and the HDFS respectively. 



liable estimate of the flux of blended sources and 3) the 
best color estimate of the sources. 

We made use of the simulated images to compute the 
completeness corrections to number counts (see §3.2) to 
systematically investigate these three different issues. The 
results of our analysis agree with those previously found 
by Martini (2000) and Saracco et al. (1999) showing that 
the SExtractor BEST magnitude (i.e. isophotal corrected 
magnitude for sources flagged by SExtractor as blended 
and Kron-like magnitude for the remaining sources) un- 
derestimates the flux when approaching the limits of the 
survey. This bias is much larger in the case of isophotal 
magnitudes (both corrected or not). We also found that a 2 
arcsec aperture photometry ('^3xFWHM in our images) 
brightened by 0.15 mag (derived by point-like sources) 
gives a good estimate of the magnitude of the faintest 
sources while underestimate the flux of Ks<20.5 sources 
and fails to recover that of blended sources. 



We thus estimated magnitudes on the basis of the fol- 
lowing steps: 

— the 2 arcsec aperture magnitude corrected by 0.15 mag 
was assigned to those "isolated" sources having an 
isophotal diameter (on the smoothed image) smaller 
than 2 arcsec. 

— the BEST magnitude was assigned to the remaining 
sources. 

This is very similar to the procedure adopted in the pre- 
vious papers (Fontana et al. 2000; Arnouts et al. 1999) 

3.2. Galaxy counts 

To derive galaxy counts, we first removed stars from the 
samples relying on the Sextractor star/galaxy classifier. 
We defined stars as those sources having a stellarity index 
larger than 0.9 and a point-source FWHM in all the bands. 
A visual inspection made on the F814W band image of the 
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HDFS confirmed the reliability of the classification. We 
found 16 and 23 stars (all of them brighter than Ks=21) 
in the CDF and in the HDFS respectively. The possibility 
that our classification fails to detect fainter stars implies 
negligible differences in the number counts at these mag- 
nitudes, since the surface density of stars is at least one 
order of magnitude lower than that of galaxies. 

We have computed the correction factors for faint, un- 
detected sources following the recipe described in Saracco 
et al. (1999) and Volonteri et al. (2000). We generated for 
each band a set of frames by dimming our final images by 
various factors and adding pure poissonian noise at the 
same level as that of the original images. Sextractor was 
then run with the same detection parameters to search 
for real sources in the dimmed frames. The correction fac- 
tor c is the mean number of dimmed galaxy which should 
enter the fainter magnitude bin over the mean number 
of detected ones. In Tab. 2 we report the raw counts n^, 
the correction factors c, the counts per square degree cor- 
rected for incompleteness and their errors a^- Errors were 
obtained by quadratically summing the Poissonian contri- 
bution = ^/inr) of raw counts, and the uncertainty 
on the correction factor Uc. 

Fig. 1 and 2 show the number counts derived for the 
Ks-band and for the J-band, respectively. We do not find 
significant deviations between the CDF and HDFS. The 
coimts follow a dlog'H / dvn relation with a slope 7j ~ 0.34 
at J> 20 and 7^ ~ 0.28 at Ks> 19. Also plotted in the fig- 
ures are number counts derived from the previous surveys 
taken from the literature. 

The surface density of galaxies we derived as well as 
the rate at which it increases with apparent magnitude, 
i.e. the slope of the counts, lie in between those derived 
previously by other surveys. 

The remarkable scatter among the surveys may be as- 
cribed to the superposition of various effects: the slightly 
different filters (K, K', Ks and J, Js), the different method 
used to estimate magnitudes, possible systematics in the 
photometric calibration and cosmic variance which could 
significantly affect the number counts at faint magnitudes 
where the area covered by the surveys is always narrower 
than few arcmin^. It is actually at magnitudes fainter than 
K=21 and J=23 that the largest discrepancies among vari- 
ous surveys occur. Finally, fields centered on high-redshift 
target objects may result in a significant excess of the 
number counts (e.g. Soifer et al. 1994). 

The CDF and HDFS Ks-band counts derived here 
arc systematically lower than those of Bersliady et al. 
(1998), Saracco et al. (1999) and McCracken et al. (2000). 
The largest discrepancy occurs with respect to Bershady's 
counts (a factor 1.7 at Ks^ 22.5) which increase faster 
(especially in the Ks band) with respect to our counts, 
the slope being ~ 0.36. Both Bershady et al. (1998) and 
McCracken et al. (2000) surveyed high Galactic latitude 
blank fields covering an area smaller than 2 arcmin^. They 
count more than 200 and 80 objects down to K~ 23.5 and 
K^ 22.5 respectively. While Bershady et al. use fixed aper- 
ture magnitudes (within ~ 2 arcsec) corrected to total on 



the basis of the object size, McCracken et al. use a Kron- 
type magnitude. Saracco et al. (1999) surveyed an area of 
20 arcmin^ centered on the NTT Deep Field (Arnouts et 
al. 1999) and including the high-redshift quasar BR 1202- 
07. They detect ~ 1000 objects to a limiting magnitude 
Kg ^ 22.5 and apply an aperture correction to the 2.5 
arcsec fixed aperture magnitude. The method we used to 
estimate magnitudes is quite similar to the one used both 
by Bershady et al. and Saracco et al. The different sur- 
face density of galaxies they found with respect to the one 
here derived may thus be ascribed to cosmic variance af- 
fecting the small area covered by Bershady et al. and the 
field surveyed by Saracco et al. due to the presence of the 
high-redshift quasar. In the case of McCracken et al. data 
their different method used to estimate magnitudes adds 
to the high count fluctuations expected in such a small 
area. 

On the other hand our counts, which agree with those 
of Moustakas et al. (1997), are systematically higher than 
those of Djorgovski et al. (1995) and of the recent very 
deep counts derived by Maihara et al. (2000) on the 
Subaru Deep Field (SDF). In particular the CDF and 
HDFS counts are in excess at magnitudes J> 23 and Ks> 
22 with respect to the SDF counts which are described by 
a slope of ^ 0.23 and they are always ~ 1.4 times higher 
than those derived by Djorgovski et al. Bershady et al. 
(1998) claim that the corrected aperture magnitude used 
by Djorgovski et al. (1995) may result in an overestimation 
of the depth of their survey up to 0.5 mag. 

It is worth noting that extrapolating the scaling rela- 
tion found for field K-band selected galaxies (Roche et 
al. 1999) to Ks~21, we found that clustering can pro- 
duce count fluctuations of ~ 25% in the magnitude range 
20.5 < Ks < 21.5 over an area comparable to the CDF 
and HDFS. These fluctuations rise up to ~ 30 — 32% over 
areas of 2-1 arcmin^. Thus, large scale structure fluctu- 
ations alone could in principle account for most of the 
observed discrepancies among the surveys at faint magni- 
tudes. 



4. Colors of galaxies 

The near-IR color of galaxies of our samples has been 
obtained by running SExtractor in the so-called double- 
image mode: J (Js) magnitudes for the Ks selected sam- 
ple have been derived using the Ks frame as a reference 
image for detection and the J (Js) image for measure- 
ments. We used the difference between the 2 arcsec di- 
ameter aperture magnitudes (~ 3xFWHM) to estimate 
the color of isolated sources while the difference between 
isophotal magnitudes has been adopted to measure the 
color of blended sources. This is to minimize the uncer- 
tainties on flux estimates occurring for blended sources 
in the case of a flxed aperture magnitude. We considered 
genuine estimates those signals exceeding 2a above the 
background while fainter signals have been replaced by 
the 2(7 upper limit. 
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Table 2. Differential number counts in the Ks-band and in the J-band derived on the CDF (upper panel) and on 
the HDFS (lower panel). Errors take into account the Poissonian error on raw counts and the uncertainties in the 
completeness correction factor c. 



In Fig. 3 and 4 the color-magnitude diagrams of the Ks 
selected CDF and HDFS samples are shown respectively: 
open circles represent galaxies detected both in the J and 
in the Ks band while the lower limits to the J-Ks color 
are represented by the arrows. It is worth noting that the 
different depth reached in the J filters (see Tab. 1) gives 
rise to a difference of ~ 0.15 mag in the limiting J-Ks color 
between the two fields. Point-like sources are marked by 
star symbols. The median J-Ks color in 0.5 magnitude bins 
for our sample (filled circle) is also plotted and compared 
to the value given in Bershady et al. (1998, filled squares) 
and Saracco et al. (1999, filled triangles). Error bars on 
the median values are the standard deviation within each 
bin. The horizontal lines represent the typical J-K color of 
a main sequence M6 star (dashed line, J-K~ 0.96) and of 
an elliptical galaxy at z > 2 (solid line). In the figures we 
also plotted the expected J-Ks color of an M* elliptical 
(solid curve), spiral (Sb, long-dashed curve) and irregular 
(short-dashed curve) galaxy in the redshift range < z < 
3. Model prediction are based on Buzzoni's (1989, 1995) 
population synthesis code: ellipticals are described by a 
single burst stellar population, spirals are modeled taking 
into account a declining star formation rate in the disk 
coupled with a spheroidal component as in the elliptical 
and irregulars follow a flat star formation rate at every 
age. 

We have compared the color distributions of the two 
samples: a feature present in the CDF sample is the bluing 



trend described by the median color occurring at magni- 
tude Ks> 18 where the J-Ks color change from ~ 1.7 at 
Ks'^ 18 to ^ 1.4 at Ks^ 21.5. The median colors agree 
well with those in Bershady et al. (1998) and Saracco et 
al. (1999) as well as the bluing trend we found, a trend 
much more evident in the case of optical-IR colors (see 
e.g. Gardner et al. 1993; McCracken et al. 2000). In the 
HDFS sample the bluing trend is less obvious. At bright 
magnitudes (Ks< 20) the CDF and HDFS samples follow 
almost the same color distribution while at fainter mag- 
nitudes the median colors in the HDFS seem to remain 
constant and arc slightly but systematically redder than 
those in the CDF. This is shown in Table 3 where we re- 
port the median colors for the two samples as a function of 
Ks magnitude and in Fig. 5 where the color distributions 
in different Ks magnitude slices for the CDF and HDFS 
samples are compared. The different behavior shown by 
the median color of the two samples is not due to the dif- 
ferent limiting J-Ks color between the two fields, at least 
down to Ks=21.5 Indeed, in this magnitude range all the 
K-band selected objects of both samples have been de- 
tected in J. 

The median color in the CDF seems to be slightly red- 
der than in the HDFS at very bright magnitudes (Ks< 19) 
while it is 0.1 bluer than the HDFS at Ks> 19. These small 
differences cannot be ascribed to systematics in photomet- 
ric calibration since most of the observations were been 
carried out during the same nights. Various effects can in 
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Fig. 3. Color-magnitude diagram for the CDF Ks selected 
sample. Open circles represent grtlaxics detected both in 
the J and in the Ks band the arrows mark the J-Ks lower 
limits. The horizontal lines represent the typical J-K color 
of a main sequence M6 star (dashed line, J-K~ 0.96) and 
of an elliptical galaxy at z > 2 (solid line) respectively. 
The expected J-Ks color for a pure galaxy population of 
ellipticals (E, solid curve), spirals (Sb, long-dashed curve) 
and irregular (Irr, short-dashed curve) types is also shown 
(Ho, qo = 50, 0.5). 
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Fig. 5. J-Ks color distributions for the whole CDF and 

HDFS sample (All) and for subsaniples selected in differ- 
ent Ks magnitude ranges. At bright magnitudes (Ks< 20) 
the two samples are described by the same color distribu- 
tion. At fainter magnitudes the HDFS sample shows colors 
slightly redder than to those of the CDF one. 

principle contribute to this difference. For instance, the 
presence of an intermediate- ^; cluster in such small areas 
would strongly influence the median colors. This could 
be the case of the CDF where the presence of a cluster 
at ^; ~ 06 — 0.7 (Cimatti, private communication) would 
affect colors at magnitudes Ks< 19, giving rise to the 
slightly redder color observed in this interval. At Ks> 19 
the larger number of red (J-Ks> 1.9) objects populating 
the HDFS with respect to the CDF (see Fig. 3 and 4) gives 
rise to the different median colors. This different number 
of red objects can be explained in terms of ERO clustering 
(see §6). Finally, the use of slightly different filters could 
produce systematics in the color distributions. The stan- 
dard J filter used to image the CDF is slightly redder and 
much wider than the Js filter (see Tab. 1) and has some 
leaks in the K band. All these features act to make the 
CDF sources bluer than the HDFS, an effect more evident 
going to high redshift where the rest-frame optical part of 
the spectra enters the J filter. We have tried to quantify 
this effect by convolving the response functions of the two 
filters with a synthetic spectrum of an Sb galaxy devoid 
of emission lines. At z = the difference between the es- 
timated magnitudes within the two filters is 0.012 mag in 
the sense that the galaxy is brighter in J (bluer in J-Ks), 
as we expected. At redshift z = 2, where the Balmer break 
enters the J filter, the difference rises to 0.083 mag. The 
possibility that an emission line falls in the J filter but not 
in the Js due to its narrower width significantly increases 
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1.44 


0.4 
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0.3 
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20.75 


1.48 


0.3 


21.25 


1.57 


0.4 


21.25 


1.40 


0.4 


21.75 


1.72 


0.4 


21.75 


1.43 


0.5 



Table 3. Median J-Ks color of galaxies as a function of 
magnitude. Errors are the standard deviation in each bin. 



this difference. For instance, the presence of a typical Ha 

emission line (which enters the J band at z ^ 0.9) would 
add 0.04 mag to the magnitude in the J filter. 

5. Comparison with models 

We have compared our number counts with different mod- 
els by using ncmod, the publicly available code by Gardner 
(1998). The aim of this comparison is not to define how 
galaxies evolve, since it is well known that number counts 
alone are not able to put severe constraints on the evolu- 
tion of galaxies. We would rather be interested in under- 
standing the most probable mechanisms able to reconcile 
optical and near-IR number counts in the light of the lat- 
est estimates of Q. 

We have considered two families of Hq = 70 Km s~^ 
Mpc~^ flat cosmological models: the first one defined by 
(f2M, f^A)=(l-0,0.0) and the second by 
(Qm, f^A)=(0.3,0.7), according to the recent results 
of the Boomerang and MAXIMA experiments (De 
Bernardis et al. 2000: Balbi et al. 2000) and in agreement 
with the results on Type la supernovae (Riess et al. 1998; 
Perlmutter et al. 1999). The expected Ks and J galaxy 
counts for the two scenarios are shown in Fig. 6 and 7 
respectively. All the models are based on the following 
assumptions: l)the local K-band luminosity function (LF) 
and mix of galaxies of Gardner et al. (1997); 2) the galaxy 
spectral energy distributions from GISSEL96 (Bruzual 
and Chariot 1993) We considered five morphological 
types of galaxies: E/SO, Sa, Sbc, Scd and Irr. We used 
in the computation the filter transmission curves of 
the VLT-ISAAC J and Ks filters. In the models, all 
galaxies formed at z = 15 with the exception of irregulars 
which are always 1 Gyr old. E/SO, Sab and Scd evolve 
following an exponential star formation rate (SFR) with 
e-folding time t = 1.0 Gyr, t = 4.0 Gyr and r = 7.0 
Gyr respectively while Scd and Irr are assumed to be 
described by a constant star formation. These are the 
basic parameters describing the pure luminosity evolution 
(PLE) "ev" models represented in the figures by the 
short-dashed curve. These models tend to underestimate 
the counts at faint magnitudes, an effect much heavier 
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Fig. 6. Galaxy count models {Hq = 70 Km s~^ Mpc~^) 
are compared to the observed Ks-band counts. The long- 
dashed and short-dashed lines represent a PLE model with 
(ev-|-du) and without (ev) absorption by dust. The addi- 
tion to the ev model of a dwarf population {a — —1.6 
for Irr) is described by the dot-dashed line (ev-|-dw) 
while the solid-line represents a model including merging 
(ev-|-dw-|-me) . The left-hand panel shows these models in 
a (f^M, f^A)=(0.3,0.7) universe, while the right-hand panel 
considers a (Om, fiA)=(l-0,0.0) cosmology. 

in the (Om, f^A)=(l-0,0.0) case. We will now explore 

the various standard possibilities to reduce the observed 
excess of faint galaxies in the near-IR with respect to our 
PLE model. 

The addition of dust extinction (long-dashed curve) 
has been introduced in modeling galaxy counts mainly 
to reduce the UV excess in galaxies and to match the U 
and B-band galaxy counts. In this case it enhances the 
observed excess at faint magnitudes even if the effect is 
minimal in the near-IR. In these models the absorption 
by dust internal to the galaxies is based on the recipe of 
Wang (1991) and Bruzual et al. (1988) who modeled the 
dust as a layer with variable thickness symmetric around 
the center of the galaxy. 

A steep faint end LF (a = —1.6) describing the latest 
types of galaxies reduces in principal the discrepancy at 
faint magnitudes. Evidence of a steep LF of field galaxies 
at infrared wavelength come from Szokoly et al. (1998) 
who find a slope a = —1.3 and from Bershady et al. 
(1999) who estimate a — —1.6. Other evidence in favor 
of a <C — 1 come from the results obtained on optically se- 
lected samples of field galaxies (Marzke et al. 1994, 1997; 
Zucca et al. 1997) and of cluster galaxies (De Propris et al. 
1998; Lobo et al. 1997; Bernstein et al. 1995; Molinari et al. 
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Fig. 7. Same as Fig. 6 but for the J-band galaxy counts. 

1998). However, sucli a steep LF, whicli contributes sub- 
stantially to the galaxy counts at optical wavelengths (see 
e.g. Volonteri et al. 2000), affects only marginally galaxy 
counts in the near-IR due to the blue optical-IR colors 
of irregular and dwarf galaxies. The small contribution to 
the counts given by this steep LF at faint IR magnitudes 
is shown in the figures by the dot-dashed curve (ev+dw), 
a model which assumes a = —1.6 for the LF of irregular 
galaxies. 

We finally considered a model including both merg- 
ing and the dwarf component (solid curve). This model 
follows the number evolution proposed by Broadhurst et 
al. (1992) who parameterized merging in the form (f)* oc 
exp[~Q/ + z)~'^ — 1)], where Q defines the merger rate 
and /3 = 1 + (2go)°-^/2 (Peacock 1987). This formalism is 
preferable to that where merging proceeds exponentially 
with (1 + z) (e.g. Rocca-Volmcrange & Guidcrdoni 1990) 
because the rate increases smoothly at high redshift and 
does not become unreasonably high at early times. 

This model fits very well the observed J and Ks galaxy 
counts both in the Qm, f^A=l-0,0.0 and Qm, f^A=0.3,0.7 
cosmology, assuming a merger rate Q=l. However, the 
same fit could be obtained without the dwarf component 
by enhancing the rate of merging to Q=1.5. 

It is worth noting that all the models without merging 
here considered predict a turnoff in the counts (at Ks^ 
20 in the qo = 0.5 models and at Ks~ 21 in the Ao = 
0.7 models), contrary to the observations. On the other 
hand, models including merging both in the qq = 0.5 and 
Ao = 0.7 cosmology provide a good fit to the data over the 
whole magnitude and wavelength ranges considered. The 
qo = 0.5 model including merging and dwarf component 
is in fact the same model previously explored by Volonteri 



et al. (2000) which provided a good fit also to the HDFS 
galaxy counts at optical wavelengths. 

The merging rate required to fit our data (from Q=l 
to Q=1.5) implies that $*(z = 1)q=i = 1.5$*(z = 0) 
and $*(z = 1)q=i.5 = 2^*{z = 0), i.e. a galaxy should 
undergo on average 0.5-1 merger events from z = 1 to 
z = 0. These rate of merging are comparable with those 
derived from the observed merger fraction based on counts 
of close pairs of galaxies and from the merger identifica- 
tion studies. The results of Roche et al. (1999) suggest 
that a galaxy will undergo from 0.4 to 0.7 merger events 
in the redshift interval < z < 1 while Le Fevre et al. 
(2000) estimate from 0.8 to 1.8 merger events in the same 
redshift range. Comparing the rate of merging derived by 
galaxy count models with that predicted by a hierarchical 
clustering scheme is less obvious. The model of Baugh et 
al. (1996) predicts that 50% of the elliptical galaxies and 
15% of the spiral galaxies have had a major merger in the 
redshift interval < z < 0.5. These numbers increase to 
~ 90% and ~ 50% at z < 1. On the other hand, a sig- 
nificant fraction of present day ellipticals (^ 40%) should 
be the result of the merging of four or more fragments 
at z ~ 0.5, while most of the spirals should contain only 
two such fragments at z ~ 1. The simple galaxy counts 
model here considered accounts for a mean merging rate 
averaged over the whole population of galaxies, which at 
z 1 should be ^ 1.5 — 2 times larger than that at z = 0. 
These numbers seem to be lower than those predicted by 
hierarchical clustering. 

In the right-hand panels of figures 6 and 7 is also shown 
for interest a PLE qo = 0.02 model (dotted-line) including 
the absorption by dust internal to the galaxies. It fits rea- 
sonably well both J and Ks-band galaxy counts. This kind 
of model, which seems to provide a good fit to the data 
over a large wavelength range (e.g. Pozzetti et al. 1996, 
1998; Metcalfe et al. 2000; McCracken et al. 2000), is in 
fact ruled out by the recent measurement of 

6. Color selections 

6.1. The fraction of z > 1 galaxies 

Kauffmann and Chariot (1998) showed that hierarchical 
galaxy formation makes substantially different predictions 
about the redshift distribution of a K-band selected sam- 
ple with respect to PLE models. In particular, hierarchical 
models predict a much smaller fraction of galaxies at z > 1 
where most of the massive galaxies should not be already 
assembled. The fraction of z > 1 galaxies in a K-band 
selected sample can thus represent a powerful test to dis- 
criminate between the different galaxy formation models. 

Fontana et al. (1999, 2000) use a photometric redshift 
technique to derive the redshift distribution of a K< 21 se- 
lected sample of 319 galaxies. They find a good agreement 
with the prediction of hierarchical models and estimate a 
fraction of 35-40% of galaxies at z > 1 and of 5% at z > 2. 

Eiscnhardt et al. (2000) suggest that the J-K color of 
galaxies is a good photometric redshift indicator, showing 
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Table 4. Lower limits to the fraction of > 1 galaxies as 
a function of Ks magnitude. 



that all but one of the galaxies redder than J-K=1.9 in 

their spectroscopic sample of ^ 50 galaxies arc at redshift 
z > 1. Such a color threshold should actually select old 
stellar systems at redshift z > 1, J-K~ 1.9 being the color 
of a present day elliptical galaxy at z > 1. Conscqiicntly, 
this color selection provides a lower limit to the fraction 
of ^; > 1 galaxies, since later- type z > 1 galaxies would 
be discarded due to their bluer color. Eisenhardt et al. 
find a lower limit of 25% for the fraction of z > 1 galax- 
ies in the EES K= 20 limited sample. This number is 
in agreement with the photometric redshift distribution 
of Fontana et al. (2000) at the same limiting magnitude. 
Indeed, a careful examination of the published photomet- 
ric redshift catalogs used by Fontana et al. (2000) shows 
that the two criteria are roughly equivalent at K< 20, as 
expected since at these bright K limits most of the selected 
galaxies are early type. 

As a first attempt to extend this test at fainter mag- 
nitudes, we have applied the selection criterion suggested 
by Eisenhardt et al. (2000) to our Ks selected samples in 
order to derive the lower limit to the fraction of z > 1 
galaxies. In the HDFS sample we applied this color se- 
lection taking into account the different .J filter used, i.e. 
restricting the selection to those sources redder than Js- 
Ks> 2 (see §4). 

In Table 4 the total number of galaxies (Ntot), the 
number of galaxies redder than J-Ks=1.9 (N2>i) and their 
fraction with respect to the total (fz>i) are shown for each 
Ks magnitude bin. We estimate a fraction of J-Ks>1.9, 
Ks< 20 objects of 14% and 7% in the CDF and HDFS 
respectively, corresponding to a surface density of 1.5 and 
0.8 objects per square arcmin. Eisenhardt et al. find a 
mean surface density of J-Ks>1.9 objects of 3.4 over the 
124 arcmin^ of the EES survey. They show that this value 
can range from 1 to 6.7 on subfields having a size compa- 
rable to the HDF. The values we find are thus consistent 
with those derived by Eisenhardt et al. and the differences 
between CDF and HDFS confirm their findings. 

At fainter magnitudes (20 < Ks < 22), the fraction 
of J-Ks> 1.9 objects rises by 13-27% in the two fields. 
These numbers are comparable to the prediction of hier- 
archical models originally computed by Kauffmann and 
Chariot (1998), and lower than the fractions predicted by 
the revised hierarchical models computed by Fontana et 
al (1999), that range from 23% of the Standard CDM to 




Fig. 8. The figure shows the J-band image (left-hand 
panel) and the Ks-band image (right-hand panel) of source 
CDF-293. Both the images are 15x15 arcsec centered on 
the source (RA=03:32:13, Dec=-27:46:41) This source is 
characterized by a color J-Ks=2.6 and an apparent mag- 
nitude Ks=20.4. 

40% of the A CDM (both at K< 21). At the same time, 
they are significantly lower than the fraction of z > 1 
predicted by PLE models, which predict up to 75% of 
galaxies in the same magnitiide range (Kauffmann and 
Chariot 1998). Taking into account the expected fluctu- 
ation due to the small size of the field, and considering 
that the .J-K> 1.9 criteria produces a lower limit to the 
actual fraction of z > 1 galaxies, we conclude that both 
the hierarchical and PLE models are still compatible with 
our findings, and that more accurate determinations of the 
redshift distributions at z > 1 are required at these faint 
magnitudes to discriminate between them. 

6.2. EROs selection at z > 2 

In §4 we derived the J-Ks color distribution of the Ks 
selected CDF and HDFS samples noting the presence 
of some sources redder than J-Ks~ 2.3 at magnitudes 
Ks> 20. This is the k-corrccted color of an elliptical galaxy 
at 2 > 1.9 or at z > 2.1, in the case of passive evolution 
(Buzzoni 1995). This color threshold would thus in prin- 
ciple select old stellar systems at z > 2 where the Balmer 
jump would be redshifted red- wards of the J filter and 
the UV emission would be depressed by the lack of young 
stellar population. On the other hand, similar colors also 
could be displayed by dusty star-forming galaxies which 
would look so red due to the heavy absorption of the UV 
emission by dust. 

In the CDF sample we count 7 sources (2% of the sam- 
ple) redder than J-Ks=2.3 at Ks< 22, i.e. a surface density 
of about 1.2 per square arcmin. In the HDFS sample we 
selected 20 sources (5% of the sample) redder than Js- 
Ks=2.4 at Ks< 22, equivalent to a surface density of 2.7 
per square arcmin. Six out of these 20 EROs lie out of 
the WFPC2 field and four on the borders. Four out of the 
remaining ten have no optical counterpart in the F814W 
band image. 

The different surface densities found are most probably 
due to the clustering properties of EROs which introduce a 
strong field-to- field variation in small areas. Indeed, Daddi 
et al. (2000) detect a strong clustering signal of EROs se- 
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Fig. 9. The J-band image (left-hand panel) and the Ks- 
band image (right-hand panel) of source HDFS-363 are 
shown. Both the images are 15x15 arcsec centered on 
the source (RA=22:32:49, Dec=-60:32:11). This source is 
characterized by a color Js-Ks=3.3 and an apparent mag- 
nitude Ks=20.8. 

lected on the basis of the R-K color from a K~ 19 limited 
sample (see also McCarthy et al. 2000). They show that 
such clustering produces fluctuations on the surface densi- 
ties of K< 19 EROs 1.7 times higher than those expected 
from pure poissonian statistic on areas of the order of ^ 25 
arcmin^. 

We selected EROs on a sample 3 magnitudes fainter 
and, in principle, on a different redshift range since we 
have used the J-Ks color selection. Thus, the results of 
Daddi et al. (2000) cannot be directly applied to our data. 
However, we could assume that the ERO clustering am- 
plitude scales with K magnitude accordingly to the am- 
plitude foimd for field K-band selected galaxies. On the 
basis of the results obtained by Roche et al. (1999) we 
thus expect that at Ks-^ 21 the ERO clustering is 1.25 
times weaker than that derived by Daddi et al. at Ks^l9. 
In the presence! of a correlation with amplitude A, the 
fluctuations a of the counts around the mean value n are 
(J — \/n(l + nAC), being C oc Area~^''^. Applying this 
relation to our data, we derive an ERO surface density 
fluctuation ct ~ 0.7 which could be consistent with the 
fluctuations we observe. 

In Fig. 8 and 9 the J and Ks-band images of the EROs 
CDF-293 and HDFS-363 are shown as an example. These 
sources have a J-Ks color of 2.6 and 3.3 respectively and 
an apparent Ks magnitude of 20.4 and 20.8. If they are 
elliptical galaxies at redshift z > 2, their k-corrected ab- 
solute magnitude would be Mks < —25, consistent with 
a passively evolved M* {z = 0) galaxy. 

Objects with extreme J-K color fainter than K> 20 
have been also found in the HDF-N by Dickinson et al. 
(2000) and in the HDF-S NICMOS field by Yahata et al. 
(2000). They discuss the possibility that they are dusty 
star-forming galaxies or ellipticals at z > 2 or z > 10 
Lyman break galaxies. Maihara et al. (2000) notice that 
most of the red sources they found in the SDF appear as 
close neighbors, thus being probable interacting systems. 
We notice that ~ 30% of the red sources we found both 
in the CDF and in the HDFS are flagged as "blended" by 
SExtractor and that the remaining 70% do not appear as 
close neighbors. A final possibility remains that some of 



the most compact sources at Ks> 21 are faint very low- 
mass stars, which can display colors redder than J-K=2 
(Chabrier et al. 2000; Dahn et al. 2000). 

The analysis of the near-IR properties combined with 
those derived by photometry at optical wavelength will be 
presented in a forthcoming paper (Saracco et al. 2001). 

6.3. Extremely blue sources 

Figures 3 and 4 reveal the presence of objects with J-Ks 
color bluer than ^0.8 mag, the color of a typical irregular 
galaxy. On the basis of the criteria we used to discriminate 
between star and galaxy (§3.2) those of them brighter than 
Ks~20 turned out to be stars. The remaining blue sources 
have been classified as galaxies since they have either a 
FWHM larger than a point source or a stellarity index 
lower than 0.9. However, almost all of them (14 in the 
CDF and 8 in the HDFS) are fainter than Ks=21 where 
star/galaxy classification is unfeasible in our data. Three 
out of the eight HDFS blue sources He in the WFPC2 field: 
while two of them are pointlike in the F814W band image, 
the remaining one displays an irregular shape. 

The first hypothesis we could make about the nature 
of these sources is that they are Galactic objects. In this 
possible explanation for the blue IR color of at least 
some of them is that they are extreme M sub-dwarfs, simi- 
lar to those studied by Leggett et al. (1998). These sources 
have metallicities as low as about one- hundredth solar, 
are characterized by effective temperatures Teff ~ 3000 
K and by masses ranging 0.09-0.15 Mq. The objects an- 
alyzed by Leggett et al. have magnitudes in the range 
9< K < 15; colors J-Ks> 0.4 mag and absolute magni- 
tudes in the range 9< Mk < 10. Most of our extremely 
blue sources are in fact in this color range. Under the hy- 
pothesis that they are M sub-dwarfs with absolute magni- 
tudes similar to those of Leggett et al., the present obser- 
vations should have reached ~ 6 kpc considering a limiting 
magnitude Ks^ 22.5. 

Sources even bluer than J-Ks=0.4 and approaching J- 
K~ can be consistent with T-type brown dwarfs (e.g. 
Nakajima et al. 1995; Cuby et al. 1999). These IR ex- 
tremely blue objects are characterized by red optical-IR 
colors. For instance, the methane brown dwarf NTTDF 
J1205-0744 (Cuby et al. 1999), the most distant reported 
to date, has a magnitude Ks=20.3. and colors J-K~ 0, I- 
K> 6. We reveal 3 faint sources in our samples with such 
a blue IR color: 2 in the CDF sample having J-Ks< 0.3 
and 1 in the HDFS sample with J-Ks< 0.1. All of them 
are fainter than Ks=22 where the uncertainty in the color 
estimate is larger than 0.3 mag. If our three sources are 
methane brown dwarfs they have to be 2.5 times more 
distant than NTTDF J1205-0744 (considering AK = 2 
magnitudes) . The analysis of the optical data will help to 
assess the nature of these sources. 

The second hypothesis about the nature of our ex- 
tremely blue sources is that they are extragalactic objects. 
Using Bruzual & Chariot models we can account for colors 
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less blue than 0.65, the case of a 1 Gyr old irregular galaxy 
at z = with a constant star formation (1 M0/yr) and a 
metallicity lower than O.2Z0. Bluer colors can be obtained 
only assuming irregulars younger than IGyr. Another pos- 
sibility is that these blue objects are z > 3 sources with 
strong UV excess, like those reported by Sullivan et al. 
(2000). They debate the possibility that the UV light 
in these objects is due to a non-thermal source such as 
QSO/AGN perhaps superimposed on a star-forming com- 
ponent. In this case we would observe the rest-frame UV 
(< 3000 A) in the J band and wavelengths short-ward of 
5500 A in the Ks band. In fact, by red-shifting the mean 
energy distribution (MED) derived by Elvis et al. (1994) 
for quasars we are not able to produce colors bluer than 
0.85 {z ~ 3). The large dispersion about this MED results 
in an uncertainty we estimate to be about 0.2 — 0.3 mag. 
However, it is worth noting that QSOs showing colors as 
blue as J-K~ 0.6 — 0.5 are observed also at redshift z < 3 
(Francis et al. 2000). 



7. Summary and conclusions 

We presented counts and colors of galaxies detected in 
the deep J and Ks images obtained with the near-IR cam- 
era ISAAC at the ESO VLT telescope centered on the 
Chandra Deep Field and on the Hubble Deep Field-South. 
The co-addition of short dithered 2.5x2.5 arcmin images 
led to a total exposure time of about 8 hours in Ks and 
to a limiting surface brightness Ks~22.8 mag/arcsec^ and 
J~24.5 mag/arcsec^ on both the fields. 

We found an excellent agreement between the counts 
derived in the two fields. On the other hand our counts 
and their slope lie in between those previously derived 
by other authors at comparable depth. In particular we 
estimate a slope of 7j=0.34 and 7/i-=0.28 at J> 20 and 
Ks> 19 respectively. We do not observe a turnoff in the 
counts down to the faintest magnitudes. 

The observed J and Ks galaxy counts can be well fitted 
in a flat universe {qo — 0.5 or Aq > 0) only assuming 
models including some degree of merging. These models 
were previously found also to give a good fit to the optical 
number counts of the HDFS. On the contrary models not 
including merging are not able to reproduce the observed 
counts, independent of the cosmology. 

Our estimated lower limits to the fraction of z > 1 
galaxies based on the J-K> 1.9 color selection criterion 
are broadly consistent with the predictions for hierarchical 
models as well as for pure luminosity evolution models. 

We identify 7 sources in the CDF (or 2% of the sample) 
and 20 source in the HDFS (or 5% of the sample) with col- 
ors redder than J-Ks=2.3 and magnitudes 20 < Ks < 22, 
equivalent to a surface density of 1.2 and 2.7 sources per 
square arcmin respectively. Their extreme red colors sug- 
gest that they are old ellipticals or dusty star forming 
galaxies at redshift z > 2 — 3. However some of the most 
compact of them could be very low-mass stars which can 
display very red near-IR colors. Another possibility, at 



least for the faintest and optically undetected ones, is that 
they are z > 10 Lyman break galaxies. 

A further analysis of their near-IR properties combined 
with those at optical wavelength will allow us to constraint 
their nature. 
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